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We study numerically the dynamics of two-dimensional vortex systems at zero temperature. In 
addition to pinned states and turbulent plastic flow, we find motion of vortices in rough channels 
along the direction of the driving force. In this decoupled channel regime we demonstrate how 
topological defects mediate the phase slip of different channels moving with different velocities. We 
thus provide important confirmation of recent analytical work describing vortex dynamics at high 
driving forces such as the moving glass theory of Giamarchi and Le Doussal. For the largest driving 
forces we find that the channels couple and observe elastic motion. 
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I. INTRODUCTION 

Vortex dynamics in the presence of disordering pinning 
show a variety of non-equilibrium ghysics and [dynamic 
phase traditions. Experiments, uQ numericalfTlIa and 
analyticaLjila work suggest that a disordered static sys- 
tem of vortices showS|-o|rdering at higher driving forces. 
Koshelev and Vinokurt3 predicted a dynamic phase tran- 
sition between plastic sliding for driving forces just above 
de-pinning and coherent motion of crystalline structures 
at high drivjing forces. Subsequently, Giamarchi and 
Le DoussalliJ predicted that the strongly driven and re- 
ordered system would be a Moving Glass, where vortices 



move elastically-coupled along static channels, such that 
they flow in the direction of the driving force along well- 
defined, nearly parallel paths in the pinning potential. 
These optimal channels (in two dimensions) or sheets (in 
three dimensions) show a roughness and are predicted 
to be a static and reproducible feature of the disorder 
configuration. _ . 

Balents, Marchetti and Radzihovsky arguedlij that in 
addition to elastically-coupled channels (no topological 
defects in the systemV-at-intermediate velocities a trans- 
verse moving smectidlSEj would exist in which motion 
of vortices in different channels is decoupled (topplpfc 
ical defects between the channels). Later workll3~t3 
mainly supported the initial findings of Giamarchi and 
Le Doussal with the addition of the moving smeptic as 
predicted by Balents, Marchetti and Radzihovsky.E_3 Dif- 
ferent names are in common use: the^ Moving Trans- 
verse GlassEj (MTG), moving smectictj and decoupled 
channelsEZI refer to the (decoupled channel motion- and 
the Moving Bragg GlassN (MBG), moving latticeN and 
coherent phaseO refer to the regime of elastically coupled 
channels. _ _ 

The theoretical descriptionsE-T 113 of these dynamic 
phases are based on elastic theory and assume either the 
absence (for the MBG) or the irrelevance (for the MTG) 
of topological defects In fact, the theory of Giamarchi 
and Le DoussaloEJ describes both regimes with the same 
equation, which is (nearly) exact for the MBG and re- 
mains an effective description for the MTG. In this work 
we investigate, for the first time, the role of topological 
defects in the MTG to check the validity of assumptions 
entering the theory of Giamarchi and Le Doussal, and 
find them to be justified. 

We review the dynamic phase diagram for a 2d vortex 
system in the presence of random disorder varying on a 
length scale much smaller than the vortex-vortex spac- 
ing. Our model describes rigid vortices in thin-films or 
decoupled pancake vortices in layered materials. We em- 
ploy a modified cut-ofO to the appropriate interaction 
force which corresponds to a logarithmic vortex-vortex 
interaction potential.Ej After annealing a vortex system 
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to zero temperature, we apply an increasing driving force 
and study the dynamics of the system systematically for 
different pinning strengths. 

Sec. fn| descri bes the simulation and the computational 
details. In Sec. Ill we present a dynamic phase diagram 
and g ive an overview of the observed dynamic phase s 
(III A). These phases are a pinned vort ex gla ss (IIIB), 
different kinds of tur bulen t plastic flow ( III C ) , a decou- 
pled c hanne l regime ([III D|) and coherently moving struc- 
tures ( HIE ). In Sec. |1V| we consider the decoupled chan- 
nel regime in detail: we report on the dependence of the 
spatial distribution of velocities in the different channels 
on the pinning landscape (IV A), we show how topologi- 



cal defects between the channels mediate the phase slip 
between channels whil e pres erving the transverse peri- 
odicity of the system (IV B), and we give information 
on the transv erse de-pinning in the decoupled channel 
regime flrVCp . Finally, we draw our conclusions in Sec. 
0. Appendix |X| contains technical information on the 
smooth cut-off used. 



II. THE SIMULATION 

A. Equation of motion 

We consider a two-dimensional vortex system and 
model the vortex motion with overdamped Langevin dy- 
namics. The total force Fj acting on vortex i is given 
by 



p vp 



Tithcrm 







(1) 



where r\ is the Bardeen-StepherO viscosity coefficient, 
Vj the velocity, F L the Lorentz force acting equally on 
all vortices, the vortex-vortex interaction, Fj p the 
vortex-pinning interaction, and F* hcrm a stochastic noise 
term to model temperature£3 The vortex-vortex interac- 
tion force for rigid vortices in thin films and pancakes in 
decoupled layers of multiJayer materials experienced by 
vortex i at position r, isEEl 
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The constant 4>o is the magnetic flux quantum, s the 
length of the vortex, /zo the vacuum permeability and 
A the London penetration depth. We employ periodic 
boundary conditions and cut off the logarithmic vortex- 
vortex repulsion potential smoothlyJij The important 
feature of this modified interaction potential is that it 
does not introduce numerical artefacts, such as topolog- 
ical defects which can result from using a naive cut-off 
potential. Details can be found in appendix ^ and Ref. 
|l9| . The cut-off distance is min(L a; /2, L y /2) where L x 
and L y are the lengths of the sides of the rectangular 
simulation cell. The lengths L x and L y are chosen such 



that a hexagonal lattice fits perfectly in the simulation 
cell. 
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A sample pinning potential. Distances are given 
in multiples of the vortex lattice spacing, oo- (a) Demon- 
stration of construction of the pinning potential in one di- 
mension: Firstly, we assign random pinning energies at dis- 
crete sites (shown as open circles) with spacing w. Secondly, 
we interpolate between those sites using cubic splines to ob- 
tain an effectively continous pinning potential. This results 
in a random pinning potential with a short-range correlator 
V(r)V(r') = g{r — r') of range w. We follow an analogous pro- 
cedure in two dimensions, (b) A part of a pinning potential 
as used in the simulations. The seven black cylinders indicate 
vortex lines separated by ao to demonstrate the length scale. 



We investigate systems with a magnetic induction of 
B = 1 T and a penetration depth of A = 1400A which 
yields a vortex density of ~ 10/ A 2 representative of typi- 
cal cuprate superconductors. The random pinning poten- 
tial we have employed varies smoothly on a length scale 
of A/25 which is of the order of the coherence length 
£. This is a representation of random pinning on the 
atomic length scale (for example due to oxygen vacan- 
cies or small clusters of oxygen vacancies) since the vor- 
tex cores effectively smooth the pinning potential over a 
length scale of the core diameter 2£. Fig. [I] (a) demon- 
strates the construction of the pinning potential in one di- 
mension. Fig. [l] (b) shows a part of the pinning structure 
used for the two dimensional system. System sizes from 
100 to 3000 vortices have been investigated. We mea- 
sure lengths in units of A = 1400A, and forces in units of 
the force fo that two vortices separated by A experience. 
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We express time in units of to = t]X/ fa-^n2TTfi X 4 /®o s 
which is in line with other simulations.E£lcJ 



B. Observables 

To distinguish different dynamic phases we moni- 
tor the topological defect density rid (defined as the 
fraction of vortices with less or more than six. near- 
est neighbors in the Delaunay triangulationtfl) and 
the distribution T(< v >) of time-averaged velocities 
<Vi> (t) — r '( t + t o)~ r '( t °) | f individual vortices i over 
time t. We also observe the structure factor of the 
system (the Fourier transform of the vortex positions), 
a measure for local hexagonal order (using bond an- 
gles 9k from the Delaunay triangulation wc compute 
= - 



^7 Efc=T d exp(i60 fc )|, where n bon d is the num- 
ber of angles in the Delaunay triangulation), the fre- 
quency spectrum of the center of mass velocity, and the 
paths of motion of vortices (two-dimensional histogram 
of vortex positions) . We create movies of time snap-shots 
of vortex positions to visualize the behavior of the sys- 
tem. 



An annealed vortex configuration is shown in Fig. 0. Af- 
ter annealing, a driving force is applied which is subse- 
quently increased every 4 • 10 4 time steps. This yields 
force- velocity characteristics which correspond to exper- 
imentally obtainable current-voltage characteristics (at 
zero-temperature). The driving force F L is related to 
the current density j via F L = sj x 3?o, and the vortex 
velocity to the induced electric field E via E = B x v, 
where B is the magnetic induction and v the vortex ve- 
locity. We investigate the modes of motion at different 
driving forces and pi nning strengths using the observ- 



ables specified in Sec. II B 



A. The phase diagram 

The different observed modes of plastic and elastic mo- 
tion are summarized in Tab. [j] (on page |5|). The second 
column of the table shows the expressions used for each 
mode of motion and a reference to the section in which 
it is described. The third and fourth columns show the 
criterion used to identify and distinguish the modes, and 
the fifth column gives further observations. 



III. THE DYNAMIC PHASES 




FIG. 2. A pinned vortex glass at pinning strength 
Ftms — l-2/o- The system has been annealed fronwnolten 
to zero-temperature. The figure shows the Delaunayf-J trian- 
gulation of vortex positions. Here, rid = 0.46, i.e. 46% of the 
576 vortices are topological defects (having more or less than 
six nearest neighbors) and are highlighted by open circles. 

Initially, we anneal the vortex system from a molten 
state to zero temperature in the presence of the random 
pinning potential. The pinning forces are obtained by 
numerically differentiating the potential. The root mean 
square value of the pinning force field is denoted by F^ s . 
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FIG. 3. Dynamic Phase Diagram of the vortex state as a 
function of disordering pinning strength, F^ s , and strength 
of the driving force (at zero temperature). 

We describe now briefly the phase diagram shown in 
Fig. ||. For weak pinning (F^ s < 0.8/o) a pinned vortex 
glass (Fig. |^) undergoes Plastic Flow (PF) and Order- 
ing P lastic Flow (OPF) for an increasing driving force 
(Sec. Ill C ). In OPF, in contrast to PF, the density of 
topological defects, rid, is lower than the density of the 
static system, n s d tatlc . We summarize PF and OPF as 
turbulent plastic flow because in both modes the motion 
of vortices is turbulent rather than laminar, i.e. the mo- 
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tion of different vortices is hardly correlated. This helps 
us to distinguish between the turbulent (chaotic) plastic 
flow of PF and OPF and the (laminar) plastic moti on of 
vortices in the decoupled channel regime (Sec. [II D , Sec. 
IV). At high driving forces the vortices move elastically 
as a coherently moving structure, i.e. e very v ortex keeps 
its nearest neighbors for all times (Sec. Ill E] ). 

For the weakest pinning strength of F^ s « JX0if 
the annealed system is a defect-free Bragg glasses and 
changes directhz-fxom the pinned Bragg glass to the mov- 
ing Bragg glassE-3 (MBG) without undergoing plastic mo- 
tion. However, since the initial configuration is annealed 
from random positions, this MBG is, in general, not 
aligned with the direction of the driving force, and the 
pinning is too weak to reorientate it. This very weak pin- 
ning regime has not been studied in detail in the frame- 
work of this investigation. 

For stronger pinning (F^ s > O.8/0) there is an inter- 
mediate regime between turbulent plastic flow and co- 
herently moving structures in which rows of vortices are 
aligned with the driving force and vortices move in pre- 
ferred channels. However, these channels are decoupled: 
vortices in different channels move with different veloci- 



ties (Sec. [V). The change from the decoupled channel 



regime to a coherently moving structure depends on the 
history of the system: in the mixed regime both modes of 
motion can be found depending on whether the driving 
force is increased or decreased. 

Related numerical work on dynamic phases has been 
performed by Moon, Scalettar and-Zimanyi,Ll Ryu et al.t3 
and Olson, Reichhardt and Nori.O It was found in Ref. 
H that as the driving force is increased, firstly the pinned 
vortex glass exhibits Plastic Flow and finally moves as 
a "moving glass" which is very likely to be the decou- 
pled channel regime. Ref. [lO] found an elastically moving 
structure with topological defects at high driving forces. 
In contrast to this work in which we have used loga- 
rithmic interactions and have varied the strength of the 
pinning forces, in Ref. [l2] the strength of an exponen- 
tially decaying vortex-vortex interaction has been var- 
ied in a system with a smaller vortex density. However, 
the results can be compared qualitatively, and Ref. [l2| 
demonstrates similar findings on Plastic Flow, decoupled 
channels and coupled channels. 



C. Turbulent Plastic Flow 

Vortices start moving if the driving force exceeds a 
critical value. We distinguish two different kinds of mo- 
tion which we refer to as Plastic Flow (PF) and Ordering 
Plastic Flow (OPF). Both types of motion show a broad 
distribution of time-averaged vortex velocities as shown 
in the inset in Fig. |4](a), which, for very small driving 
forces, has another peak at zero velocity. We call the 
motion OPF if the density of topological defects, rid, is 
below the defect density, n d tatIC , the system would have 
if no driving force was applied. Otherwise we call it PF 
(Tab. §. 

We observe PF for driving forces just above the critical 
de-pinning force. The topological defect density is higher 
than for the static system because some vortices are sta- 
tionary and others are squeezing past them. Within the 
PF regime we find two modes of motion: For driving 
forces just above the de-pinning current we find a bi- 
modal distribution in the time-averaged vortex velocity 
showing a peak at zero velocity. Thus, there are some 
vortices that are permanently pinned (at least over the 
simulated time). By contrast, for higher driving forces, 
whilst at any one time some vortices may be stationary, 
no vortices are permanently pinned. These data con- 
firm earlier findings of Spencer and JensenEil employing 
a simpler model. For clarity, Fig. ^ does not distinguish 
between these two types of PF. 

In the OPF regime, where the topological defect den- 
sity, rid, is lower than for the static system, we observe 
that the instantaneous velocity distribution shows no 
peak at zero velocity (i.e. all vortices are in motion), 
whereas this is not the case for both types of PF de- 
scribed iiij-the last paragraph. Faleski, Marchetti and 
MiddletonEj used the term crinkle flow to describe mo- 
tion of vortices via correlated displacements of patches 
of vortices. Our observations of OPF suggest that the 
definition used here for OPF (rid < n d tatlc ) is equivalent 
to the definition of crinkle, flow introduced by Faleski, 
Marchetti and MiddletoncJ (absence of a peak at zero 
velocity in the instantaneous velocity distribution). 



D. Decoupled channels 



B. Pinned vortex system 

For sufficiently small driving forces the system is 
pinned and the velocity distribution shows a single peak 
at zero velocity. For pinning strengths above w 0.04/o 
we see the number of topological defects increasing with 
pinning strength and no long range order exists. We thus 
refer to the pinned system as a vortex glass, and such a 
configuration is shown in Fig. ^. 



For sufficiently strong pinning and intermediate driv- 
ing forces (Fig. ^) we find that vortices arrange in lines 
orientated along the direction of the driving force (Fig. 
[|). These lines move with different velocities in the direc- 
tion of the driving force. This type of motion is described 
in detail in Sec. IV and is called decoupled channel mo- 
tion. 
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Name 


Criterion 


Observations 




rid 


Plastic 
modes 


Plastic Flow (PF) and some vor- 
tices permanently pinned Sec. 1IIC 


Broad, and 
peak at zero 


n d > n| tatlc 


Turbulent flow, system partly pinned 


Plastic Flow (PF) and no vortices 
permanently pinned Sec. [II C 


Broad 




Turbulent flow, peak at zero in instantaneous ve- 
locity distribution (i.e. some stationary vortices) 


Ordering Plastic Fiow (OPF) 
Sec. QIC 


Broad 


, ^static 

n d < n d 


Turbulent flow, no vortices have zero velocity 
in instantaneous velocity distribution "crinkle 
motion" {i.e. all vortices moving) 


Decoupled channels, 

Moving Transverse Glass (MTG) 

oec. |L VI 


Separated 
5-peaks 


r\ static 

< rid <C n d 


Motion in uncoupled channels in direction of driv- 
ing force, topological defects between channels, 
critical transverse force 


Elastic 
modes 


Coherently moving structure 
without defects. Moving Braee 
Glass (MBG), Sec. |[IIe| 


Single 5-peak 


n d = 


Motion in coupled channels in direction of driv- 
ing force, washboard frequency in noise spectrum, 
critical transverse force 


Coherentlv moving structure with 
defects, Sec. HIE 


Single 5-peak 


0<7id<n d 


Vortices generally aligned with the direction of the 
driving force 



TABLE I. Overview of observed plastic and elastic m odes of motion. r(<u>) is the distribution of time averaged vortex 
velocities, and rid is the density of topological defects (Sec. |ll Ej ) . The topological defect density of the annealed system without 
any applied driving force is n d tatlc . 
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FIG. 4. The decoupled channel regime. Periodic boundary conditions are applied in x and y directions, (a) Velocity 
histogram of time-averaged individual vortex velocities F(<v>) for a driving force F£ = 5.0/o- Inset: same for turbulent 
plastic flow with F% = 0.9/o. (b) Delaunay triangulation of one time step with the time-averaged individual vortex velocity 
plotted in the third dimension. There are four distinct groups labeled A,B,C and D of vortex-channels travelling along the 
^-direction with different velocities, (c) Two-dimensional view of (b) and topological defects are highlighted, (d) Change 
between initial (•) and final (x) vortex positions in the frame of reference of one of the vortices in group C. 
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x x 

FIG. 5. The Moving Bragg Glass. Left: Histogram of vortex positions. The driving force is acting from left to right along 
the a;-direction, and vortices move in rough channels, like beads on a string. The inset shows a slightly enlarged version of the 
channels and positions of vortices for one time step are shown as circles. Right: Delaunay configuration of one snap shot of the 
same system. Although the channels in the left plot are rough, there are no topological defects in the Moving Bragg Glass. 





FIG. 6. A coherently moving structure with 6 dislocations (12 topological defects). 
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E. Coherently moving structure 

We observe two different kinds of coherently moving 
structures: (i) either a Moving Bragg GlassliJ (MBG) 
which is free of topological defects (Fig. ||), or (ii) a 
hexagonal system similarly aligned with the direction 
of the driving force but with a few dislocations (Fig. 
||). Both configurations move elastically, i.e. each vor- 
tex keeps its nearest neighbors for all times. 

For strong pinning (F^ a > O.8/0) and increasing 
driving force the transition from the decoupled chan- 
nel regime to a coherently moving structure results in 
a MBG if the groups of cou pled c hannels have the same 
vortex line density (see Sec. IVB). If the groups of cou- 
pled channels have different line densities then the dis- 
locations between them are frozen into the coherently 
moving structure. For weak pinning (F^ s < O.8/0) the 
vortices do not move in decoupled channels for intermedi- 
ate driving forces, and the system changes directly from 
Ordering Plastic Flow (OPF) to a coherently moving 
structure. Again, elastically moving systems with and 
without topological defects are observed. Our data from 
simulating current-voltage characteristics with increasing 
driving force suggest that the configuration at high driv- 
ing forces is usually a hexagonal system aligned with the 
driving force with a few dislocations. However, a MBG 
configuration is occasionally observed. 

The coherently moving structures we observe are al- 
ways aligned with the direction of the driving force for 
pinning strengths > O.8/0. For smaller pinning strengths, 
configurations develop occasionally which are ncj 
with the driving force. It has been arguec 
that this alignment minimizes power dissipation, and 
our results axe. in agreement with other numerical 
investigations!]!!!] in which the high-velocity configura- 
tions are generally aligned with the driving force. 

For a MBG we find peaks at multiples of the washboard 
frequency luq = 2ir <v cm > /oq in the Fourier spectrum of 
the center-of-mass velocity v cm (t) of the system, where 
<> denotes a time- average and ao is the lattice constant 
of the vortex lattice. Whereas this temporal periodicity 
is not existent for the velocity of an individual vortex, we 
also find it in the energy of the system. Clear peaks in 
the Fourier spectrum can be observed up to frequencies of 
w lOOwo- For a single particle the washboard frequency is 
observable if it slides through a periodic potential. Here, 
we have a random potential, but a periodic system. One 
thus finds the washboard frequency in observables that 
depend on all vortices, such as the center-of-mass veloc- 
ity or the energy, but not for individual vortices. The 
wasfabpard frequency has been found experimentally in 
AC2a~Eil and DCE3 measurements, and numericallytHl in 
a similar regime. 

In summary, we observe occasionally a MBG at high 
driving forces, but most of the final configurations are 
hexagonal systems aligned with the driving force with 
a few dislocations pairs. However, it could well be that 



finite temperatures or larger systems would favor the cre- 
ation of a MBG at high driving forces: as yet, it is not 
clear what is the "dynamic ground state" of these sys- 
tems. Our data cannot be used to decide whether a 
MBG exists in two-dimensions or whether the MTG is 
the only stable phase,t3'E3 since for a system of a given 
size if the velocity is sufficiently large then all channels 
couple and appear to be a MBG. Another open question 
is whether periodic boundary conditions can favor a re- 
ordering of a disordered vortex system.LJ The exploration 
of these questions is computationally expensive though 
new methods for evaluating interactions in the system 
may make this feasible.ESEJ 



IV. DECOUPLED CHANNELS 

This section describes a plastic mode of motion which, 
due to its quite different properties, is separated from 
the Sec. Ill C on turbulent plastic flow. As visible in the 
dynamic phase diagram in Fig. y the decoupled channels 



are only observed for sufficiently strong pinning. 

Increasing the driving force (which acts in the re- 
direction) from the turbulent plastic flow regime further, 
transforms the time-averaged velocity distribution from a 
broad peak as observed for turbulent plastic flow to sev- 
eral clearly distinct peaks as shown in Fig. |] (a). In ^ (b) 
four different velocity levels are visible, each of these cor- 
responding to one peak in the velocity histogram. Thus, 
vortices move in four groups of coupled channels and, 
within a group, all channels travel with a constant ve- 
locity in the direction of the driving force. Plot |J (c) 
shows the lattice structure of (b) in a two-dimensional 
projection. Vortices with more or less than six nearest 
neighbors are highlighted by a grey shade. We see that 
the groups of coupled channels are separated from each 
other by one 5-7 dislocation (a pair of vortices one having 
5 and the other 7 nearest neighbors). 

Plot H (d) shows the initial (•) and final ( x ) positions 
of vortices in the frame of reference of one of the vor- 
tices in group C, and the initial and final positions are 
connected by a straight line, demonstrating that vortices 
never change the channels in which they move. This is 
a particularly interesting point since the moving glass 
theoryO assumes that the topological defects (which are 
hard to treat analytically) between groups of coupled 
channel do not destroy transverse periodicity. Thus, the 
observation that these topological defects do not intro- 
duce chaotic motion of vortices and that the defects just 
decouple the different channels supports the theory of 
Giamarchi and Le DoussalJia 

A series of runs shows that generally for larger driving 
forces the number of channels (and thus the size of each 
group) which are coupled and move with the same veloc- 
ity increases, until the system shows elastic motion. For 
decreasing driving forces the number of coupled channels 
decreases until each group of coupled channels exists of 
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only one or two channels. For even smaller driving force, 
the systems exhibit turbulent plastic flow. 

For the mixed regime shown in Fig. || we find that the 
mode of motion depends on the history of the system: In- 
creasing the driving force for a system in the decoupled 
channel regime into the mixed regime results in motion 
in decoupled channels. On the other hand, reducing the 
driving force for a Moving Bragg Glass (MBG) into the 
mixed regime, yields elastic motion. For driving forces 
above the dotted line in Fig. || all systems show elastic 
motion and below the dash-dotted line all systems show 
smectic motion. The data suggest that both the MBG 
and the decoupled channels are metastable steady states 
which are separated by an energy barrier. In future we 
will explore whether finite temperatures are able to over- 
come this barrier. 




FIG. 7. This plot shows for the decoupled channel regime 
in which areas of the simulation vortices prefer to travel. On 
a 700x700 cell grid, a two-dimensional histogram of vortex 
positions has been created. The darker a cell, the more vor- 
tices have been counted within that cell over the duration of 
the simulation. 

Fig. ^ shows an accumulation of vortex positions using 
a grid of 700 x 700 cells. It demonstrates that the channels 
in which vortices move are not strictly static but slightly 
broadened (see, for example, j/k3), although vo rtices 
never change channels. Further analysis in Sec. JlVB 
shows that the 5-7 dislocations highlighted in Fig. |J (c) 
move with time in the x-direction parallel to the driving 
force. Presumably this requires slight corrections of the 
static channels, which results in their blurred form visi- 
ble in Fig. 0. This is supported by results given in Fig. || 
which show that for the MBG in the absence of disloca- 
tions the resulting channels are strictly static (and not 
blurred) . This may indicate that the theoretical modeH 
predicting strictly static channels for the MTG may be 
too simple. 
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FIG. 8. Intensity of two Bragg peaks for the Moving 
Transverse Glass shown in Fig. ^. The (01) peak is relative 
constant at a value of approximately 0.85 and it measures the 
order along the y-direction. In contrast, the peak (10) mea- 
sures also order along the :r-direction, and oscillates strongly. 
The variations are due to the different groups of coupled chan- 
nels sliding past each other. The deviation of the (01) peak 
from 1.0 (as for a perfect lattice) is due to the roughness of 
the channels. The inset shows the lattice vectors used to label 
the peaks. 

Fig. H shows the square modulus of the structure factor 
for two k- vectors for the MTG. The large value of the 
(Ol)-peak indicates the transverse order of the system. 
The small and oscillating (10) peak shows the strongly 
reduced order in the x-direction. A time Fourier trans- 
form of the (10) signal reveals the frequencies with which 
the different groups overtake each other by one lattice 
spacing. This is exactly what is expected for a MTG in 
a finite system, and is in agreement with the theoretical 
predictiont3~t2l for the smectic regime that any peaks in 
the structure factor with a non-zero x-component should 
vanish in an infinitely large system. 

The data shown are in agre emen t with the theoreti- 
cally predicted moving smectictSll^l which is also called 
Moving Transverse Glasslia (MTG), and decoupled chan- 
nel regimeO A MTG has previously been identified by 
a numerical study of Olson, Reichhardt and Norijij and 
the Delaunay triangulation of a snap-shot of their system 
looks qualitatively like Fig. [| (c). Kolton, Dommgup. 
and Gr0nbech- Jensen numerically found smectic states,E3 
and earlier the .-numerical studies of Moon, Scalettar 
and T. ZimanyiD on moving vortex systems suggested 
the possibility of phase slips of different channels. Fur- 
ther new results on the MTG are presented in the next 
two sections concerning the spatial distribution of vortex 
chan nel ve locities depending on the pinning landscape 
(Sec. IV A), the mechanis m of uncoupled channels slid- 
ing past each othe r (Sec. IV B ), and the transverse de- 
pinning (Sec. IV C] ). 
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A. Dependence of the spatial velocity distribution 
on pinning landscape 

We find a correlation between the particular pin- 
ning potential employed (representing details of the mi- 
crostructure causing vortex pinning in the material) and 
the positions and velocities of the different groups of cou- 
pled channels. It should be noted that, although for the 
data shown in Fig. ^ (b) the vortices with the maximum 
velocity are located in the central region of the sample 
(y ~ ± 1.5), this is not an edge effect: for other samples 
the maximum is located at different y-positions. Remark- 
ably, in both the decoupled channel regime (Fig. ^) and 
the turbulent plastic flow regime (not shown here) the 
fastest flow is located in the same part of the simulated 
material. ■— ■ 

It has been suggested!^ that the different velocities ob- 
served in the decoupled channel regime and shown in 
Fig. H may be rolated to the experimentally observed 
fingerprint effect .□ In fact, it seems that for both the 
turbulent plastic flow and the decoupled channel regime 
the same areas of the pinning potential allow for bet- 
ter (or worse) pinning. This is not obvious since in the 
turbulent plastic flow regime vortices flow more or less 
individually along highly tortuous paths whereas in the 
decoupled channel regime they move in a much more cor- 
related way. 



B. Channel sliding mechanisms 

Fig. |J shows that the number of vortices per line of 
vortices (the line density) differs from group to group 
by 1, and exactly one dislocation between the groups is 
required to accommodate-this difference ([| c). In Fig. ^ 
a time series of DelaunayEfl triangulations of snap-shots 
of a part of Fig. |] is shown, demonstrating how a moving 
dislocation allows group C to move faster than group D. 

Fig. [l^ shows this process for two neighboring chan- 
nels having the same longitudinal vortex density. Again 
the phase slip is mediated by dislocations which travel 
along the channel. However, since initially there are no 
dislocations between the channels a dislocation pair is 
created. These dislocations travel away from each other, 
allowing the phase slip between the upper and lower two 
lines. When the dislocations meet again (due to periodic 
boundary conditions) they annihilate. 

Only the two mechanisms shown in Fig. [)] and [lO] have 
been observed. The situation with the same longitudinal 
vortex density has been observed less frequently which 
may indicate that this is energetically more expensive. 
However, in macroscopic systems the two mechanisms 
described are less distinguishable and may coexist: the 
local vortex line density differs around each of the dis- 
locations in Fig. 10. Thus, the process shown in Fig. O 




may just be a more detailed study of the phase slipping 
process in Fig. [HI for each of the dislocations. 



FIG. 9. Dislocation mediated phase slip in decoupled 
channel regime for different longitudinal vortex densities 
showing snap shots of 9 time steps. The open circles in each 
plot mark two fixed vortices, and the upper one defines the 
frame of reference. The Lorentz force acts from left to right. 
The upper two lines of vortices visible in each plot belong to 
group C in Fig. ^ and the lower to lines to group D. Thus, 
in the relative frame of reference, the lower two lines move 
to the left. The two black filled circles indicate topological 
defects in each snap-shot and represent a dislocation. These 
mediate the phase slip while moving to the left as can be seen 
by comparing the open circles in plot 1 and 9. 



FrjOm a figure in the work of Olson, Reichhardt and 
NoritL3 we identify varying longitudinal vortex densities 
and Burgers vectorscH parallel to the driving force, both 
in agreement with our results. 

In conclusion, the detailed mechanism of decoupled 
channels moving past each other has been identified for 
the first time. The phase slip is mediated by disclina- 
tion pairs which either exist between separate groups of 
coupled channels with locally different line densities of 
vortices, or the disclination pairs are created dynami- 
cally and in pairs when sufficient shear stress has built 
up. These results may help in finding a starting point 
for a theoretical description of the dynamics of dislo- 
cations, such as a Kosterlitz-ThoulessEa theory for non- 
equilibrium systems. 
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plot 12 (time s1ep=271 000) 



FIG. 10. Dislocation mediated phase slip in decoupled channel regime for equal longitudinal vortex densities. The open 
circles mark four fixed vortices and black filled circles show topological defects. Here, in the relative frame of reference, the 
lower two lines move to the right. The phase slip is realized by two dislocations consisting of a 5-7 and a 7-5-disclination pair 
(counting first the nearest neighbors for the upper disclination) . These emerge from a pair-anti-pair creation process in plot 2. 
The 5-7-dislocation on the left moves to the left and the 7-5-dislocation to the right-hand side. In this process they allow the 
lower lines to pass a lattice spacing to the right. Between plot 6 and 7 the 7-5-dislocation leaves the simulation cell at the 
right-hand side and enters it again on the left-hand side. Finally, in plot 10, the two dislocations meet again and annihilate 
each other. The topological defect in the lowest row is not important here. 
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C. Transverse de-pinning 

Following the theory of Giamarchi and Le Doussal,0 
at zero temperature, a non-analytical response of the vor- 
tex system to a small transverse force is expected for the 
moving glass, i.e. for the MBG and the MTG. In agree- 
ment with our results in both the MBG and the MTG the 
existence of such transverse barriers have been observed 
in simulations BE3E3'e3 The transverse de-pinning of the 
MBG has recently been describedE3 and here we report 
on the transverse de-pinning of the MTG. We have found 
that the transverse de-pinning of a MTG may happen in 
two ways: 

Mechanism A: At a certain strength of the transverse 
force (below the transverse de-pinning force), some vor- 
tices change the rows in which they have moved so far, 
such that after this first step the different groups of cou- 
pled channels A, B, C and D (in Fig. I|) have the same 
number of vortices in each row. Then, all these rows move 
with the same velocity, i.e. the system has changed from 
a MTG to a MBG. This MBG does not (yet) move in 
the transverse direction. Only when the transverse driv- 
ing force is increased even further, the system de-pins 
in the transverse direction, and moves elastically in the 
longitudinal and the transverse direction as described in 
Ref. H 

Mechanism B: At the transverse de-pinning force the 
system rearranges plastically such that the rows of vor- 
tices (which are aligned with the rc-axis in Fig. ^) become 
orientated with an angle to the x-axis after the change. 
The new direction of the rows is not the same as the di- 
rection of the total driving force (adding the small trans- 
verse force to the main driving force along the x-axis) . 

We have found that the critical transverse force is 
higher for mechanism A. Our early investigations have 
shown that the transverse de-pinning of the MTG is an 
intricate matter and further studies are required to reveal 
under which circumstances mechanism A or B appears. 



V. CONCLUSIONS 

We have modelled the dynamics of vortices in two 
dimensions using overdamped Langevin dynamics with 
a logarithmic vortex-vortex—interaction potential which 
includes a modified cut-oftJ to avoid introducing nu- 
merical artefacts into the simulation. We have com- 
puted a dynamic phase diagram as a function of pinning 
strength and driving force. We find pinned vortex sys- 
tems, different kinds of turbulent plastic flow, and for 
large driving forces motion of vortices in rough chan- 
nels along the direction of the driving force. Depend- 
ing on pinning strength and driving force the motion 
in different channels can either be coupled or decou- 
pled. These phases can be identified with the predicted 
Moving Bragg GlassM (MBG) and the Moving Trans- 
verse GlassEj (MTG) as described in recent theoretical 



models.0-0 

We have studied the MTG in detail and report on the 
dependence of the vortex channel velocities on the pin- 
ning landscape. We have identified how topological de- 
fects mediate the phase slip between channels moving 
with different velocities, and we have shown that vor- 
tices never change the channels in which they are mov- 
ing, i.e. the dislocations in the system do not produce 
chaotic motion of vortices, thus preserving transverse pe- 
riodicity. Together with the observed critical transverse 
force for the MTG and the MBG in these simulations, 
our findings strongly support_the moving glass theory 
of Giamarchi and Le Doussallia which assumes that the 
dislocations in the MTG do not introduce additional ef- 
fects that may destroy transverse periodicity (and thus 
the critical transverse force) in the MTG. 

Our findings may also, help in finding an extension to 
the Kosterlitz-ThoulcssH theory for non-equilibrium sys- 
tems. 
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APPENDIX A: SMOOTH CUT-OFF 

We employ a smooth cut-off for the vortex-vortex in- 
teraction following the ideas described in Ref. Here 
we give details on the particular interpolating function 
we have chosen. 

The vortex-vortex interaction has to be cut off for dis- 
tances greater than a cut-off distance, b. Assume the 
interaction force is given by /(r). For short-ranged in- 
teractions it is sufficient to use an interaction /(r) which 
is /(r) for r < b and zero otherwise: 



f(r) 



m 





r < b 
r > b. 



For long-ranged forces this approach results in artificial 
configurations .113 However, those problems can be over- 
come by reducing f(r) smoothly to zero near the cut-off 
distance b. One needs to introduce another distance, a, 
and a polynomial p(r), such that a < b and that p(r) 
interpolates between /(a) at a and and zero at b: 



f(r) 




r < a 

a < r < b 

r > b. 



It is requiredllj that f(r) shows C 1 continuity at a and b, 
and its derivative at b to be zero: 



dr 
dp 
dr 



f(a) =p(a), 
p(b) = 0, 
dp 
dr 



0. 



(Ala) 
(Alb) 

(Ale) 
(Aid) 



r—b 



n 



We have used a third order polynomial 



p(x) = ax 1 — C3X 3 + C2X 2 + c\x + co 



i=0 



and the coefficients Cj are completely determined by equations (Ala) to (Aid). Writing f'(r) for one finds 




(a - 6) 3 



V 



f'(a)a-f(a)b-2f(a) \ 
-f'(a)a 2 + 3 f(a)(a + b)- af'(a)b + 2 f'(a)b 2 
(2 f(a)a 2 - a/'(a)6 - 6 f(a)a - f'(a)b 2 ) b 
-f'(a)b 2 a 2 + 3 f(a)ab 2 + f'(a)ab 3 - b 3 f(a) J 



The cut-off distance b is determined by geometrical constraints (see Sec. II A). We follow Ref. |19| and choose the 
distance b — a over which the interaction is reduced to zero to be three lattice spacings, so that a — b — 3ao. In this 
work f(r) oc 1/r. Fig. 11 shows a schematic plot of the smooth cut-off and the interpolating polynomial. 

To compute the potential energy of the system it is required to integrate — p(x) to represent the smoothed interaction 
potential for a < r < b. The integration constant is determined by requiring continuity of the interaction potential at 
r = a. 
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FIG. 11. Demonstrating the shape of the interpolating polynomial p(r) (thick line) which smoothly reduces the vortex-vortex 
interaction force f(r) to zero. For clarity we have chosen f(r) = 1/r. The interpolation starts at the fading distance a = 7.5oo 
and reduces the interaction force to zero at the cut-off distance b — 10.5ao, where ao is the average vortex lattice spacing. See 
text for details. 




12 



1 R. Thorel, J. Phys. (Paris) 34, 447 (1973). 

2 S. Bhattacharya and M. J. Higgins, Phys. Rev. Lett. 70, 
2617 (1993). 

3 U. Yaron, P. L. Gammel, D. A. Huse, R. N. Kleiman, 

C. S. Oglesby, E. Bucher, B. Batlogg, D. J. Bishop, K. 
Mortensen, K. Clausen, C. A. Bolle, and F. De La Cruz, 
Phys. Rev. Lett. 73, 2748 (1994). 

4 M. J. Higgins and S. Bhattacharya, Physica C 257, 232 
(1996). 

5 M. C. Hellerqvist, D. Ephorn, W. R. White, M. R. Beasley, 
and A. Kapitulnik, Phys. Rev. Lett. 76, 4022 (1996). 

6 M. Marchevsky, J. Aarts, P. H. Kes, and M. V. Indenbom, 
Phys. Rev. Lett. 78, 531 (1997). 

7 F. Pardo, F. de la Cruz, P. L. Gammel, E. Bucher, and 

D. J. Bishop, Nature 396, 348 (1998). 

8 A.-C. Shi and A. J. Berlinsky, Phys. Rev. Lett. 67, 1926 
(1991). 

9 K. Moon, R. T. Scalettar, and G. T. Zimanyi, Phys. Rev. 
Lett. 77, 2778 (1996). 

10 S. Ryu, M. Hellerqvist, S. Doniach, A. Kapitulnik, and D. 
Stroud, Phys. Rev. Lett. 77, 5114 (1996). 

11 S. Spencer and H. J. Jensen, Phys. Rev. B 55, 8473 (1997). 

12 C. J. Olson, C. Reichhardt, and F. Nori, Phys. Rev. Lett. 
81, 3757 (1998). 

13 A. E. Koshelev and V. M. Vinokur, Phys. Rev. Lett. 73, 
3580 (1994). 

14 T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 76, 3408 
(1996). 

L. Balents, M. C. Marchetti, and L. Radzihovsky, Phys. 
Rev. Lett. 78, 751 (1997). 

16 L. Balents, M. C. Marchetti, and L. Radzihovsky, Phys. 
Rev. B 57, 7705 (1998). 

17 S. Scheidl and V. M. Vinokur, Phys. Rev. E 57, 2574 
(1998). 

18 P. Le Doussal and T. Giamarchi, Phys. Rev. B 57, 11356 
(1998). 

19 H. Fangohr, A. Price, S. Cox, P. A. J. de Groot, G. J. 
Daniell, and K. S. Thomas, J. Comput. Phys 162, 372 
(2000). 

20 J. R. Clem, Phys. Rev. B 43, 7837 (1991). 

21 J. Bardeen and M. J. Stephen, Phys. Rev. 140, A1197 
(1965). 

22 P. M. Chaikin and T. C. Lubensky, Principles of condensed 
matter physics (Cambridge University Press, Cambridge, 
1995). 

23 A. van Otterlo, R. T. Scalettar and G. T. Zimanyi, Phys. 
Rev. Lett. 81, 1497 (1998). 

24 M. de Berg, M. van Kreveld, M. Overmars, and O. 
Schwarzkopf, Computational Geometry (Springer Verlag, 
Berlin, 1997). 

25 T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 72, 1530 
(1994). 

26 M. C. Faleski, M. C. Marchetti, and A. A. Middleton, Phys. 
Rev. B 54, 12427 (1996). 

27 A. Schmid and W. Hauger, J. Low. Temp. Phys. 11, 667 



(1973). 

28 J. Miillers and A. Schmid, Ann. Phys. (Leipzig) 4, 757 
(1995). 

29 S. Bhattacharya, J. P. Stokes, and M. J. Higgins, Phys. 
Rev. B 43, 1835 (1991). 

30 J. M. Harris, N. P. Ong, R. Gagnon, and L. Taillefer, Phys. 
Rev. Lett. 74, 3684 (1995). 

31 J. M. Harris, K. Krishana, N. P. Ong, R. Gagnon, and L. 
Taillefer, J. of Low Temp. Physics 105, 877 (1996). 

32 Y. Togawa, R. Abi ru, K. Iwaya, H. Kitano, and A. Maeda, 
cond-mat/0009476 . 

33 V. Vinokur, (private communication). 

34 S. J. Cox, G. J. Daniell, H. Fangohr, and A. M. Robinson, 
submitted. 

35 A. B. Kolton, D. Dommguez, and N. Gr0nbech- Jensen, 
Phys. Rev. Lett. 83, 3061 (1999); Phys. Rev. B 62, R14657 
(2000); Phys. Rev. Lett. 86, 4112 (2001). 

36 S. Bhattacharya, (private communication). 

37 N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(Saunders College, Orlando, 1976). 

38 J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 
(1973). 

39 C. J. Olson and C. Reichhardt, Phys. Rev. B 61, R3811 
(2000). 

40 H. Fangohr, P. A. J. de Groot, and S. J. Cox, Phys. Rev. 
B 63, 064501 (2001). 



13 



